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Abstract 



Recent results concerning the internal structure of static spherically- 
■ symmetric non-Abelian black holes in the Einstein- Yang-Mills (EYM) 

. theory and its generalizations including scalar fields are reviewed and 

discussed with an emphasis on the problem of a generic singularity 
in black holes. It is argued that in the theories admitting a violation 
of the naive no-hair conjecture the structure of singularity is essen- 
(D ', tially affected by the "hair roots" . This invalidates an image of a 

1^ ' non-Abelian black hole as a Schwarzschild black hole sitting inside 

the soliton. We give an analytic description of the generic oscillatory 
^ , approach to the singularity in the pure EYM theory in terms of a 

^ I divergent discrete sequence and show that the mass function is expo- 

nentially growing "in average" . The second type of a generic approach 
to the singularity in hairy black holes is a "power-law mass inflation" 
which is realized in the theories including scalar fields. Both singular- 
ities are spacelike and no Cauchy horizons are met in the full interior 
region in conformity with the Strong Cosmic Censorship conjecture. 
Black holes violating this conjecture exist only for certain discrete val- 
ues of the event horizon radius thus forming a subset of zero measure. 
PACS: 04.20.Jb, 97.60.Lf, ll.15.Kc 
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1 Introduction 



The no-hair conjecture has played an important role in understanding 
the nature of black holes. As a result of the early investigation of various 
field theories coupled to gravity it was generally believed that the only fields 
which may extend outside the event horizon (apart from external ones not 
related to the black hole itself) are those associated with the conserved or 
topological charges carried by the hole. 

Later a notable counter example to the naive no-hair conjecture was found 
in the framework of the Einstein- Yang-Mills theory. Although no classical 
glueballs may exist in the flat spacetime because of purely repulsive nature of 
the constituent vector fields, such particle-like objects become possible when 
gravitational attraction is taken into account (Bartnik-McKinnon (BK) so- 
lutions 0). It was shown that the lower mass BK particle topologically 
is similar to the sphaleron of the Weinberg-Salam theory, with a substan- 
tial difference, however, due to existence in the BK case of the gravitational 
negative mode. Historically, just the gravitational instability of the BK so- 
lutions was discovered first |^ and later invoked as an argument favoring 
the sphaleron interpretation of the BK solutions 0. More detailed analysis 
showed, however, that the gravitational (even parity) negative modes of the 
BK particles had nothing to do with an expected sphaleron picture. Genuine 
sphaleronic features of BK solutions are manifest via the odd-parity negative 
modes along the vacuum-to- vacuum direction in the functional space 
Physically the lower mass BK solution may be seen as interpolating between 
the neighboring topologically distinct YM vacua in the EYM theory. It car- 
ries neither electric nor magnetic charge, while the Chern-Simons number 
1/2 is naturally associated with it 0, 01- 

It was soon realized that the EYM theory also admits the black hole 
counterparts to the BK solutions 0, ^. The SU{2) EYM black holes form a 
two-parameter family labeled by a continuously varying radius of the event 
horizon r/j and the number n of oscillations of the YM field in the strip [—1, 1] 
bounded by the neighboring YM vacua values. Since no asymptotic charges 
are present, these black holes apparently violate the no-hair principle. Other 
solutions sharing the same property were found in the non-Abelian gauge 



theories including scalar fields [10, 11, 12] as well as in the gravity coupled 



Skyrme model. Meanwhile, a common feeling reasonably persisted that this 
kind of hair should be distinguished from the "allowed" hair associated with 
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the global charges. It was suggested [|T3[ to interpret a non-Abelian field 



structure outside the event horizon as a "wig" rather than a genuine hair. 
Indeed, the EYM black hole can lose its YM hair as a result of the gravita- 
tional instability leaving a bare Schwarzschild black hole. As far as the black 
holes inside the magnetic monopoles are concerned, they were often thought 
of (at least for the horizon radius much smaller than the monopole size) as 
tiny Schwarzschild black holes sitting inside the solitons [Q . If so, the inside 
singularity should remain unaffected by the wig. 

However, little was known until recently about the realistic internal struc- 
ture of non-Abelian black holes (apart from some qualitative discussion in 
) . Closer investigation has shown that the idea of a Schwarzschild black 
hole sitting inside regular solitons is essentially incorrect. The interior struc- 
ture and the character of the singularity inside the EYM black hole with an 
arbitrary (continuously varying) radius of the horizon are strikingly different 
from those described by the Schwarzschild geometry. Both the Schwarzschild 
and Reissner-Nordstrom type singularities (predicted in |^) are encountered 
indeed, but only for some discrete values of the horizon radius. This "sec- 
ond quantization" in non-Abelian black holes follows mathematically from 
the same kind of a non-linear boundary value problem (now in the interior 
region) as the "first quantization" of the BK particle mass. In this sense the 
picture of a Schwarzschild black hole inside the BK particle covers only the 
set of solutions of zero measure in the parameter space. Although it could be 
anticipated by continuity that the non-Abelian hair should penetrate inside 
the event horizon, an essential role of hair in the formation of the singular- 
ity seems not to have been clearly realized before. A generic singularity in 
the spherically symmetric t-independent EYM coupled system turns out to 
be spacelike and of the oscillatory nature as was often suggested in view of 
the Bianchi IX analysis by Belinskii, Khalatnikov and Lifshitz (BKL) 
In the black hole context, however, the cosmological counterpart is given 
by the (non-Bianchi) closed Kantowski-Sachs cosmology, and the nature of 
oscillations observed is rather different from that in the BKL case. 

Our results [|I^ concerning the EYM generic black hole interior solutions 



were confirmed by Breitenlohner et al. [T^ who also attempted to extend 
the analysis to monopole black holes in the EYM-Higgs (EYMH) theory 
with the triplet Higgs. In the latter case no oscillations were observed nu- 
merically, and the behavior of the mass function near the singularity was 
found to be monotonous. This behavior was interpreted in |]T^ as exhibit- 
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ing an exponential "mass inflation" . Soon after we have found an analytical 
solution for the mass function [|18| in the SU{2) EYMH theories both with 
the (real) triplet and (complex) doublet Higgs using a consistent truncation 
of the system of equations near the singularity. Such a truncation is possible 
also in the EYM-dilaton (EYMD) theory in both cases the behavior of 
the metric near the singularity is dominated by the kinetic (gradient) con- 
tribution of the scalar field. In terms of the radial dependence of the mass 
function one finds a power-law divergence towards the singularity. Therefore 
no mass-inflation (in the usual sense pop is detected inside the EYMH black 
holes where the singularity is dominated by the scalar field. A similar scalar 
dominated asymptotic behavior was found earlier in the framework of the 
Kantowski-Sachs cosmology pT| driven by a non-linear scalar field (in the 



Abelian model). In this latter case no black hole counterpart may exist be- 
cause of the no-hair theorem for a non-linear one-component scalar field, but 
locally the solution is the same. Our formula for the power-law mass function 
in the scalar-dominated regime was reproduced later in |^ (see this volume) 
though without necessary restrictions on the power index following from the 
consistency of the truncation. It is worth noting that the domains of variation 
of the power index (depending on the horizon radius of the black hole) are 
different in the EYMH and EYMD cases, so in the Kantowski-Sachs interpre- 
tation the corresponding singularities look differently. The scalar-dominated 
singularity of the "power-law mass-inflationary" type inside spherical black 
holes in various field models including scalar fields seems to be a rather gen- 
eral phenomenon. This singularity is spacelike in conformity with the Strong 
Cosmic Censorship. Both oscillatory and the power-law singularities in hairy 
black holes therefore support this principle by the genericity argument. 

The plan of the talk is as follows. First we discuss the local solutions near 
the singularity which may be seen as dressed Schwarzschild and Reissner- 
Nordstrom singularities. Then (Sec. 3) we give an analytic description of 
the oscillatory approach to the singularity in the generic EYM black holes. 
In Sec. 4 it is shown that no exponential mass inflation develops inside the 
EYMH and EYMD black holes where the typical power-law divergence of 
the mass function is met. We conclude with a brief discussion of a new role 
of the black hole hair which was revealed in the investigation of black hole 
interiors. 
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2 Hairy Schwarzschild and 

Reissner— Nordstrom singularities 

We start with the pure SU{2) EYM system 

S=^l{-R + F'}V^d'x, (2.1) 

where F is the SU{2) field, assuming the metric to be invariant under the 
time translations and the SO (3) spatial rotations: 

A 2 2 

ds^ = ^de - ^dr^ - r\de^ + sin^ ed<p'') . (2.2) 

Here the metric functions A, a depend only on r. This parameterization of 
spacetime is suitable until a turning point of r is reached. It happens that for 
asymptotically fiat solutions there are no such points neither in exterior, nor 
in interior regions up to the curvature singularity at r = 0. Such points exist, 
however, for asymptotically non-fiat solutions, in which case another chart 



should be used . Here we will not be dealing with a priori asymptotically 
non-fiat solutions, so the coordinates ( |2.2| ) are good both in the exterior 
and interior regions. Of course one should keep in mind that in the region 
where A < the radial coordinate r is timelike, while t is spacelike, so 
that t-independence means spatial homogeneity of the spacetime rather than 
staticity. 

As usual, we choose the t-independent spherically symmetric magnetic 
ansatz for the YM potential 

A={W{r)-l){T^de~Te?>med^), (2.3) 

where T^p^Q are spherical projections of the SU{2) generators. The value 
W = 1 is a trivial YM vacuum, while W = —1 corresponds to the vacuum 
of the next topological sector (this state is related to the trivial vacuum by 
a "large" gauge transformation). 

The field equations consist of a coupled system for W, A 

Af/'+ 1--W'=--, (2.4) 



-y + 2At/2 = l-^, (2.5) 
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where V = (W^ — 1), U = W'/r, and a decoupled equation for a: 

-^lna = U\ (2.6) 

These equations admit black hole solutions in the domain r > rh for any 
radius of the event horizon |@] . The solutions for W outside the horizon lie 
within the strip — 1 < < 1 approaching ±1 asymptotically. They are spec- 
ified by the number n of nodes of W thus forming a discrete set for each r^,. 
Local solutions in the vicinity of the regular event horizon of a given radius 
contain one free parameter W{rh) 0. A "quantization" of W{rh) results from 
an imposition of the boundary conditions at infinity. Each asymptotically 
fiat exterior solution starts with some uniquely determined value W^"(r/i) at 
the horizon. All other W{rh) 7^ W^{rh) generate asymptotically non-fiat so- 
lutions. This quantization does not lead to the quantization of the black hole 
mass since the radius of the horizon still remains an arbitrary continuously 
varying parameter. 

To start the analysis of the internal structure of non-Abelian black holes, 
which is an essentially numerical problem, one has to explore first local so- 
lutions of the field equations near the singularity. Two such local branches 
were found earlier [0. The first one is the Schwarzschild-like (S) solution, it 
corresponds to the vacuum value of the YM field |Vr(0)| = 1. Introducing 
the mass function m(r), A = — 2mr, one has 

W = -1 + br^ + - 86)rV(30mo) + O(r^) , 

m = mo(l - + 86V^) + 26V^ + O(r^) , (2.7) 

where mo, b are (the only) free parameters. This local branch is non-generic 
already by counting free parameters (a generic solution of the system ( p.4|, 
p.5| ) should have three free parameters). 

The second is the Reissner-Nordstrom type branch which can be found 
assuming the leading term of A to be a positive constant (related to the 
charge parameter = A(0)). This requires W{0) = Wq 7^ ±1 and gives 

W = Wo- Wor^ {2Vo) + cr^ {2Vo) + O(r^) , 

A = V^- 2mor + + 2Woic + moWo/V^y + O(r^) , (2.8) 

what corresponds to the RN metric of the mass mo and the (magnetic) charge 
= Vq, Vq = V^(VFo)- The expansion contains three free parameters Wq, 
mo, c (i.e., is locally generic). 
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These two local solutions may be regarded as describing "hairy" Schwarz- 
schild and Reissner-Nordstrom singularities. Their essential distinction from 
the usual Schwarzschild and Reissner-Nordstrom singularities consists in 
presence of the additional free parameters {b in the first case and c in the 
second) responsible for hair degrees of freedom. 



We have also found the third local power series solution |]T5| assuming a 
negative value for A(0) (i.e., imaginary P): 

W = Wo±r- lVorV(2K)) + O(r^) , 

A = T m^V^r + 0(r2) , (2.9) 

o = ai(r2 T ^WoryVo) + ©(r^) . 

Here there is only one free parameter {Wo) for W, A. The corresponding 
space-time near the singularity is conformal to i?^ x S*^: after a time rescaling 
one obtains 

ds^ = r\dr^ - de - de^ - sin^ Od^^) . (2.10) 

This geometry was encountered in the previous study of black hole interiors in 
the framework of the perturbed Einstein-Maxwell theory ^ and called 
homogeneous mass- inflation model (HMI). 

It is easy to realize that neither of these asymptotics may correspond to 
a generic black hole. Imposing boundary conditions in the singularity, we 
obtain the same kind of the singular boundary value problem as one encoun- 
tered in the exterior problem where a similar role is played by the asymptotic 
flatness condition. This interior boundary value problem leads to the second 
quantization condition, now for the event horizon radius r^. Therefore, the 
EYM black holes with the S and RN type interiors may constitute only the 
zero measure set in the whole EYM black hole solution space. 

The system (|2.4H2.5|) was integrated numerically in the region < r < 
using an adaptive step size Runge-Kutta method for various Vh = 10~*, 10^. 
The integration started at the left vicinity of the event horizon where the 
local power series solution contains one free parameter Wh = W{rh) satis- 
fying the inequalities \Wh\ < 1 and 1 — < Vh which are the necessary 
conditions for the asymptotic flatness [§. For given Vh, the interior solutions 
meeting the expansions ( p.7| - |2.9| ) may exist only for some discrete Wh- A 



numerical strategy used to find such W{rh) consisted in detecting the change 
of sign of the derivative W'. In the S-case we found the curve W{rh) which 
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Table 1: S- and RN-type solutions. 





S-type, n = \ 


RN-type, n = 2 


RN-type, n = 3 




0.613861419 


1.273791 


1.0318420 




-0.8478649145 


-0.113763994 


-0.10185163 


^_ 




0.02171654 


0.08948446 


W{<d) 


-1 


-1.212296124 


-1.3566052 


a(0) 


0.2263801 


5.991210 X 10-3 


1.751928 X 10^3 


Mass 


0.8807931 


1.018002 


1.000277 



starts at — 1 as r/j ^ and approaches —0.1424125 for large (Fig. |l|) (with- 
out loss of generality we choose Wh < 0). Our S-curve intersects the n = 1 
branch of the family of trajectories W^{rh) corresponding to the set of ex- 
ternal asymptotically flat solutions. Parameters of this black hole are shown 
in Tab. |l|, its global behavior is depicted in Fig. ^ (for higher n S-solutions 
see [0). 

Interior solutions of the RN-type, meeting the expansions ( |2.8|) in the 
singularity, were found for rh > = 0.990288617. The corresponding curve 
W{rh) (also shown in Fig. 0) intersects the trajectories W^{rh) for all n > 
2. These solutions possess an inner Cauchy horizon at some r_ < Vh with 
|W^(r_)| >1 (Fig.i. 

Solutions of the third type ( p.9|) were studied numerically starting from 
the vicinity of the origin. The unique solution has been found for the 
horizon data subject to the necessary conditions for an asymptotic flatness 
\Wh\ <l, l-Wl<rh for the upper sign in and W{Q) = -0.9330656, 
corresponding to Vh = 1.889088. This solution, however, does not meet any 
value W^"(r/j) and thus does not represent a black hole. Thus a pure HMI 
interior can not be attributed to the EYM black holes. But it turns out to 
be an unstable flxed point of an asymptotic (truncated) dynamical system 



1^. Deviations from this fixed point correspond to generic EYM interior 



solutions sharing the same property to possess no Cauchy horizons. 

More general families of internal solutions not restricted by the asymp- 
totical flatness condition were found in |jl^ (in particular, other internal 
solutions of the third type). They do not have, however, a direct signiflcance 
for the EYM black holes. 
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Figure 1: W{rh) for the S- and RN-type interior solutions. Dashed hues — 
W'"{rh) for n = 1,2 (higher-n curves he between the n — 2 one and the 
boundary Vh — I — Wh^, dotted hne). Note that S and RN curves W{rh) do 
not merge. 
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Figure 2: The n = 1 EYM black hole (S-type). EH — events horizon. 
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3 Oscillatory approach to singularity 



Since both hairy Schwarzschild and hairy Reissner-Nordstrom singularities 
are encountered only for certain discrete values of rh (and hence the black 
hole masses), while the external solutions exist for continuously varying rh-, 
we must look for an alternative regime of approach to the singularity for 
generic r^- It is generically observed during numerical integration from the 
events horizon towards the origin that sooner or later (depending on and 
Wh) the metric function A starts to oscillate in the negative region with a 
very fast growing amplitude |jl5| . Because of huge numbers encountered in 



these oscillations the accuracy of the computation can hardly be maintained, 
so one has to seek appropriate truncations of the system of equations to 
describe the asymptotic regime analytically. Numerically it is observed that, 
when oscillation progress, the right hand side of the Eq. ( p.4| ) becomes small 
with respect to terms at the left hand side. Neglecting it one obtains the 
following approximate first integral of the system: 

Z = A[/a/r = const, (3.1) 

which relates oscillations of the mass function to the evolution of a. Numer- 
ical experiments also show that while the YM function W remains almost 
constant up to r = 0, its derivative is still non-zero and is rapidly changing on 
some very small intervals of r. The function U exhibits a step-like behavior 
being constant with high accuracy during almost all the chosen oscillation 
cycle (Fig. ^ and then jumping to a greater absolute value corresponding to 
the next cycle. It is clear from ( |2.6|) that a is exponentially falling down with 
decreasing r while U ~ const, whereas in the tiny intervals of [/-jumps a 
remains almost unchanged. So a tends to zero through an infinite sequence 
of exponential falls with increasing powers in the exponentially decreasing 
intervals. Combining this with ( p.l|) and the above mentioned properties of 
U one can deduce rather detailed description of the metric behavior. 

Let us denote by the value of radial coordinate where A has fc-th 
local maximum. Soon after passing this point, the function U stabilizes at 
some value Uk approximately equal to the doubled value at the point of local 
maximum (similarly, U increases by about a factor of two when approaching 
the local maximum, whereas A is almost stationary). Then, according to 
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0.2 0.4 0.6 0.8 1 1/10 



Figure 4: The first oscillation cycle for an EYM solution. The functions A 
for EYM (solid line) and EYMD (dashed) are shown in lower half-plane. In 
upper half-plane — mass functions m(r) (analogously) and the function U for 
EYM (dotted). All functions are power rescaled with the power index 1/10. 
Here rn = A; Wh = -0.283993 for EYM, Wh = -0.298357, (j)h = 0.05623 for 
EYMD (asymptotically flat solutions with one node of W.) 
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( |2.6| ), a is equal to 



Ulir^ - rl 



(3.2) 



From (P.ll) one finds that, while f/^ ~ const, 



A(r) = r exp \ul(rl - r^)] . (3.3) 

This function falls down with decreasing r until it reaches a local minimum 

at 



1 ^, v/|A(r,)| 

In what follows, in view of the observed fact that in the course of oscillations 
of A the YM function W changes insignificantly, we will put V = const. 

Therefore, the mass function is inflating exponentially while r decreases 
from rk to Rk- After passing Rk, an exponential in ( |3.3| ) becomes of the 
order of unity, hence A starts to grow linearly, and the mass function m(r) 
stabilizes at the value = m{Rk). Such a behavior holds until the point 
of local maximum of A/r^ is reached; this takes place when A ^ —V"^ at the 
point 

r 2 

prfecxp -{UkTk) ■ (3.5) 



' |A(r.)r 

After this a rapid fall of |A| is observed causing a violent rise of \U\. Then 
the term 2AU'^ in the Eq. ( p.5| ) becomes negligible and consequently at this 

stage 

UA ^ -V^Uk , (3.6) 

while r practically stops. This implies that very soon A reaches the next local 
maximum at the point r^+i ^ r^, while m{r) rapidly falls down to rrik+i- At 
the point of local maximum of A one has in the Eq. ( |2.(j| ) | A| ^ V'^, then in 
view of the smallness of r we find 

mn)\^ I ^, ■ (3.7) 



To obtain the estimates by the order of magnitude we will neglect all 
numerical coefficients elsewhere except for the power indices in exponentials. 
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in particular putting U{rk) = Uk, and omitting also (quasi-constant) factors 
V. With this accuracy one obtains from (|3.3|)-(|3r^): 



r,+i = M,-^ , rl^, = RkRu+i , = ^ exp , (3.8) 



|A(.)H(^) 



f Rk\'^ fk+1 



, , ^ exp 

Tk J r-k Rf 



2 



Thus, introducing a variable Xk = [rk/RkY (^ 1)? we can derive the following 
recurrent formula 

Xk+i = xfe^' , (3.10) 

which shows that is an exponentially diverging sequence. In terms of Xk 
one has 

^-^ = Xke-^^/\ (3.11) 
rk 

this relation can also be understood as a ratio of the neighboring oscillation 
periods since ^ Tk+i- Values of the function |A| at the points rapidly 
tend to zero: 

|A(r,)|=x,-\ (3.12) 

so we deal with an infinite sequence of "almost" Cauchy horizons as r ^ 0. 
At the same time the values of |A| at the points Rk grow rapidly: 

\A{Rk)\ = xf/'e^-/\ (3.13) 
and the values of the mass function grow correspondingly as 

^-x,-V^/^ (3.14) 



Ml 



k-l 

While r decreases form to Rk, the function a rapidly falls down to the 
value (Tfc = a{Rk), which then remains unchanged until the point r^+i. As 
the singularity is approached, the sequence ak decreases according to 

^ = e-^^/2_ (3.15) 

Therefore for a generic (continuously varying) r/j one observes a rather unex- 
pected approach to the singularity through an infinite sequence of oscillations 
with an exponentially growing amplitude. 
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Apart from the discrete picture given above one can find a truncated two- 
dimensional dynamical system clearly showing an infinitely oscillating nature 
of the generic solution. As we have already noted, the oscillation region 
starts with an exponential fall of A which typically occurs after passing 
a local maximum r™"^ (Fig- §)> and the right hand side of ( p.4|) becomes 
comparatively small with respect to other terms. Another important feature 
is that the YM function W (contrary to A) possesses a finite limit Wq in the 
singularity. Omitting the right hand side of (|2.4|) , replacing W by its limiting 
value Wo, and neglecting 1 as compared with V'^/r'^ one can derive from the 
system ( p.4| - p.5| ) the following two-dimensional dynamical system 



q = p, 

p= (3e-«-l)j9 + 2e"2'?-l/2, (3.16) 

where A = — (Vq^/2) exp(g), and a dot stands for derivatives with respect to 
T = 21n(r/j/r). This system has one (focal) fixed point {p = 0, q = In 2), 
corresponding to some imaginary charge RN-like local solution of the type 
( |2.9|) , with eigenvalues A = (l±zVT5)/4. Its phase portrait is shown in Fig. ^ 
together with an invariant set p = —e''^ — 1/2 corresponding to the RN-type 
local solution (|2.8| ). The generic oscillating solutions lie above this curve. 
The phase motion in this region is unbounded, and there are no limiting 
circles. One can easily show that the rotation in the phase plane never stops 
and the limit q = — oo (A = 0) can not be reached. The metric function A 
remains negative valued as r — > and passes an infinite sequence of local 
maxima and minima. The class of oscillating solutions is stable in a sense 
that a small deviation from one such solution moves us to another oscillating 
solution. Thus the existence of the consistent two-dimensional truncation of 
the full set of equations proves that the oscillating approach to the singularity 
is generic for the EYM black holes. The dynamical system ( p.l6|) was derived 



in our paper [113 1 and reproduced in a slightly different notation in [ff^] 



4 Power- law "mass- inflationary" singularity 

In non-Abelian theories including scalar fields, such as EYMH and EYMD, 
another regime of the generic approach to the singularity is realized. In 
this case the scalar kinetic term becomes dominant soon after entering an 
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Figure 5: The beginning of A-oscillations for n = 1 EYM BH solution, 
rh = 2.4, W{rh) = -0.31652531. EH — events horizon. 
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Figure 6: Phase portrait of the dynamical system (3.16), RN — an invariant 
set, corresponding to the RN-type solution (dashed — zero slope lines). 
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asymptotic regime what results in a monotonous behavior of the metric. 
Consider first the SU{2) EYMH theory 

S = ^ J !^~R-F' + 2\D<!>\'-^m'-vr]v^9d'x, (4.1) 

where $ is a Higgs field in either vector (real triplet) or fundamental (complex 
doublet) representations, D$ is the corresponding YM covariant derivative 
(in the doublet case |D<l>p = (D^Y^D^), |<l>p = $^<l>) and without loss of 
generality both the Planck mass and the gauge coupling constant are set 
to unity. In the flat space-time the triplet version of the theory gives rise 
to regular magnetic monopoles, the doublet version — to sphalerons. New 
physically interesting configurations emerge when gravity is coupled in a self- 
consistent way, in particular, static spherical black holes exist in both cases: 
monopole 0] and sphaleron ||l2l . 



Static spherically symmetric configurations of the YM fields are still given 
by the ansatz ( p.3|) , while the Higgs field is $"Ta = 0(r) Tr in the triplet case, 
and $ = (j){r)v in the doublet one, where v is some (here irrelevant) spinor 
depending only on the angle variables. In both cases 0(r) is the only real 
scalar function of the radial variable. 

The system of equations following from ( [4. 1|) with this ansatz may be 
presented as a set of three coupled equations for W, 0, and A = — 2mr: 

r'^ / r 1 oW r 

^) +A0'2 = l-2V-g, (4.3) 

(A0')' + Ar0'3 = (4.4) 
where Q = 2AW''^/r'^, and 

V ^ V{W, 0, r) = ^ + ^(0^ - vr + P\ (4.5) 

with P = W(j) in the triplet case and P = {W + in the doublet one. An 
equation for a now becomes 

(Ina)' = -W'^ + r(j)'^, (4.6) 
r 
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and can be easily integrated once W, are found. 

Local expansions of three types listed in the Sec. 2 can be easily gener- 
alized to include Higgs IJI^, and the same arguments can be used to prove 



that neither of the corresponding global solutions is generic. Meanwhile, the 



forth local branch can be found analytically in this case |T^. It can be de- 
rived using the consistent truncation of the non-linear systems of equations. 
An appropriate truncation consists in omitting from the equations all matter 
terms except for those related to the gradient of the scalar field. This reduces 
to dropping the right hand sides in Eqs. (|4^ - |4.4|) . The resulting system can 



be easily disentangled leading to the following decoupled equations: 

W 

W" - — = 0, (4.7) 

r 

0" + - = , (4.8) 

r 

which can be solved as follows: 

W = Wo + br\ (4.9) 

(f) = (f)Q + klnr . (4.10) 

Here Wq, (pQ, b, k are free parameters. Note that contrary to the expansions 
discussed in the Sec. 2, this is not a power series solution. Higgs field is 
logarithmically divergent, so that its derivative diverges as r~^, this is why 
the corresponding terms become dominant for sufficiently small r. Once W 
and (f) are found, the metric function A can be obtained by integrating the 
simple equation 

A' 1 

- = --r0'^ (4.11) 
A r 

what gives 

A = -2mor(^-^') (4.12) 

with the fifth (positive) constant mo. Hence, by counting free parameters, 
this is a generic solution with non-positive A. Now, to find whether our 
truncation, one has to substitute the solution into the Eqs. (|4.2| -|4l^) and to 
check whether the right hand side terms are indeed comparatively small. One 
finds the following condition: k'^ > 1. This means that the metric function 
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A is divergent at the singularity. The corresponding mass-function is also 
divergent according to the power-law 

m=^. (4.13) 

The asymptotic behavior of a dominated by the scalar term then reads 

a = aor^' , (4.14) 

with (positive) constant ctq. 

This local solution can be interpreted as exhibiting a "power-law mass 
inflation" . As a matter of fact this regime has nothing to do with the usual 
(exponential) mass-inflation which presumably takes place once a Cauchy 
horizon is approached. 

The second example of the scalar-dominated singularity is given by the 
EYMD theory: 

S = ^j{-R + 2(V0)2 - e-^^F^} ^gd'x , (4.15) 
The equations of motion for W, A, now take the form 

AU' - 2AU(j)' = WV/r - J^W , (4.16) 
(A/r)' + A0'2 = ^-2Af/V2<^, (4.17) 
(A0')' + Ar0'3 = J^ - 2A(0V + l)f/2e-2'^ - 1 , (4.18) 

where T = 1 — y^e~^'^r~^ , V = W"^ — 1 . The remaining equation for a 
reads 

(lna)' = r(0'2 + 2f/V2'^) . (4.19) 

Similarly to the EYMH case for sufficiently small r the right sides of the 
Eqs. ( |4.1(j| )-( |4.18| ) become small in comparison with the left hand side terms 
and one gets the following truncated system 

(lnf/)'-20' = 0, 

[ln(A/r)]' = [ln(A0')]' = • (4.20) 

Its integration gives the following five-parameter (i.e., generic) family of so- 
lutions 

W = Wo + 6r2(^-^) , 

A = -2/ir(^-^') , (4.21) 
d) = c + \n (r^^ 
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with constant Wq, b, c, /i, A. The vahdity of the truncated equations ( [4.20| ) 
now can be checked by substituting the asymptotic solution ( [4.21|) into the 
full system ( |4.16D -( P:.18| ). For consistency it is sufficient that the following 
inequahties hold: 

V2-1<A<1, (4.22) 



which is in agreement with the numerical data. 

From ( [4.21| ) it follows that the mass function diverges as r 
to the power law: 

ji_ 



The corresponding a tends to zero as 



air 



air 



according 
(4.23) 

(4.24) 



where ai = const. A typical EYMD solution is shown in Fig. ^. 

The only difference with the EYMH case is that the region of the power 
index ( [4. 221 ) is different. This leads to different picture of the singularity in 
the Kantowski-Sachs interpretation: a point like singularity in the EYMH 
case, but a cigar singularity in the EYMD case 



5 Discussion 

Non-Abelian black holes turned out to be a useful laboratory to explore the 
nature of the singularity inside "realistic" black holes. Unlike the gravity cou- 
pled Abelian field models such as Maxwell, Maxwell-dilaton-axion, or more 
general string-inspired systems, non-Abelian models have an advantage to 
give rise to several qualitatively different possibilities as far as the singular- 
ity is concerned. Moreover, in spite of the absence of exact analytic solutions, 
one can find a correspondence between the singularity structure and external 
parameters such as the black hole mass. Thus the relative weight of different 
type interior structures in the parameter space may be found. This opens a 
new way to probe the Strong Cosmic Censorship using the genericity argu- 
ment. An advantage of this approach is that a non-trivial information may 
be extracted already at the level of static ("eternal") black hole solutions. 
One can speculate that non-Abelian electroweak theory should replace the 
Maxwell electrodynamics at high energies which are reached in the course of 
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the mass-inflation inside a (perturbed) Reissner-Nordstrom black hole due to 
the phase transition similar to that in cosmology. Therefore a further fate of 
the black hole should be in the scope of a non-Abelian theory. Although the 
full time-dependent picture is likely to be quite complicated, it is reassuring 
to feel that the static spherical non-Abelian black holes choose the spacelike 
singularity in conformity with the Strong Cosmic Censorship. 

The new features observed in the interiors of non-Abelian black holes 
are due to the field components which are not connected with the conserved 
charges and are usually termed as a hair. Violation of the no-hair conjec- 
ture in black holes traditionally was attributed to their external appearance. 
New results clearly show that the hair is equally important for the internal 
structure of black holes. Hair "roots" penetrate up to the singularity sup- 
plying it with additional degrees of freedom. A tiny black hole inside a large 
magnetic monopole generically has a very different internal structure than 
the Schwarzschild black hole. Hairy black holes have hairy singularities. 

Qualitatively speaking, the interiors of the EYM black holes look like the 
hair-perturbed Schwarzschild or Reissner-Nordstrom interiors. So it is not 
surprising that one encounters an exponential growth of mass whenever the 
metric reaches an "almost" Cauchy horizon. This phenomenon, first noted 
ISH and termed as "mass inflation" in [jl^, is, however, only a half of 
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the story. In the oscillating regime one observes two qualitatively different 
"mass-inflations" : the first is the local inflation at some very short interval of 
each oscillation cycle, the second is associated with the exponential growth 
of mass from cycle to cycle |19] while the singularity is approached. 

If scalar fields are present, the dynamics gets a new dimension. Once 
the scalar component is excited, it soon becomes a dominant factor of the 
evolution which changes drastically the approach to the singularity. In the 
scalar-dominated regime no mass- inflation is manifest, instead one observes 
a behavior which we call a "power-law mass- inflation" , i.e., a power-fashion 
divergence of the mass-function near the singularity. This behavior has a 
different origin as compared with the usual mass-inflation, in particular, it is 
not related to an approach to the Cauchy horizon. Rather it follows from the 
coupled Einstein-scalar field dynamics which is essentially Abelian. Indeed, 
the same behavior near the singularity was earlier observed in the Kantowski- 
Sachs cosmology with a non-linear (one-component) scalar field source. Still 
the non-Abelian nature of the model is substantial, otherwise the (asymp- 
totically flat) static black holes are prohibited by the no-hair theorems. 
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The analysis given here is purely classical, a few words are in order about 
the relevance of the results to the full quantum theory. Vacuum polariza- 
tion of the conformal scalar field on the HMI background was considered in 



2J]. It was found that the correction to the mass function diverges more 
strongly than in the classical case. Hence the singularity is not smoothened 
but rather intensified. In the oscillating regime there is no hope to compute 
quantum effects quasi-classically. Moreover, huge values of the mass function 
(in Planck's units) encountered soon after entering such a regime indicate 
that the quantum behavior of the model should be considered nonpertur- 
batively and may well be qualitatively different from the classical picture. 
However, the conclusion about the spacelike nature of the generic singularity 
is unlikely to be changed. 

D.V.G. thanks the Organizing Committee for hospitality and support 
during the Workshop. The research was supported in part by the RFBR 
grants 96-02-18899, 18126. 
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